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In a previous paper we have considered the Regge-Wheeler equation for fields of spin 
s = 0, 1 or 2 on the Schwarzschild spacetime in coordinates that are regular at the 
horizon. In particular, we have constructed in Eddington-Finkelstein coordinates exact 
solutions in terms of series that are regular at the horizon and converge on the entire 
open domain from the central singularity to infinity. Here we extend this earlier work 
in two different directions. Firstly, we consider in Eddington-Finkelstein coordinates a 
massive scalar field that can serve as a dark matter candidate. Secondly, we extend the 
treatment of the massless case to Painleve-Gullstrand coordinates, which are associated 
with radially infalling observers. 

PACS numbers: 04.70.Bw, 03.65.Pm, 04.30.Nk, 11.80.Et 


1. Introduction 

Under the assumption that the back-reaction to the spacetime geometry can be ne¬ 
glected, the propagation of scalar, electromagnetic and gravitational perturbations 
is described by linear wave equations that have to be solved on a given background 
spacetime. For this article the background is chosen to be the Schwarzschild space- 
time. The investigation of wave propagation in such a black-hole spacetime gives 
insight into typical phenomena like interference patterns and scattering, (quasi-) 
normal modes and black hole evaporation. In a recent article 1 we have shown how 
coordinate systems that are regular at the horizon can be used for constructing ana¬ 
lytic solutions to the (massless) Regge-Wheeler equation. In particular, we have used 
Eddington-Finkelstein coordinates for constructing such solutions in terms of series 
that are convergent on the entire open interval r £ ]0, oo[ . We will use these results 
here and extend the treatment by i) considering additionally a massive scalar field, 
which brings a new parameter into the game, and ii) constructing solutions for the 
perturbation equations in the observer-associated Painleve-Gullstrand coordinates. 

Considering a massive scalar field establishes a connection to cosmology since 
it may serve as a dark-matter model.^ We have chosen to present this treatment 
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in Eddington-Finkelstein coordinates, because it is the best adapted coordinate 
system for this situation and allows to penetrate the horizon 1 following the ingoing 
radiation. Giving the solutions to the general case of a massless bosonic field with 
spin s on the other hand allows to investigate black hole physics in coordinates of a 
radially infalling observer and is therefore, at least in principle, connected to a real 
physical measurement situation. 

The article is organized as follows: We start with some preparations in section 
[2] where, after introducing our unit conventions, we recall the basic features of two 
different coordinate systems for the Schwarzschild spacetime, namely Eddington- 
Finkelstein and Painleve-Gullstrand coordinates. Furthermore, the confluent Heun 
equation that will be the central differential equation in this work is introduced and 
some of its properties are briefly discussed. 

In section [3] we consider a massive scalar field on a Schwarzschild background 
spacetime and give the wave equation that describes the propagation of this field in 
Eddington-Finkelstein coordinates. We reduce this wave equation to a purely radial 
differential equation and show how this radial equation can be transformed to a 
confluent Heun equation. Thereupon, we give the generalized Riemann scheme and 
the local solutions at the horizon in terms of the standard confluent Heun function 
HeunC. 

In section |4] we consider a general massless bosonic spin s held in Painleve- 
Gullstrand coordinates that are associated with observers which follow radial time¬ 
like geodesics. Starting from the wave equation in this coordinate system we arrive 
at a purely radial equation that belongs to the class of confluent Heun equations 
as well. We construct local solutions at the horizon in terms of the confluent Heun 
function and asymptotic series at spatial infinity. These local solutions have conver¬ 
gent representations only in a bounded domain of convergence, but we examine an 
analytic continuation procedure to overcome this problem. 


2. Preparation 

2.1. Eddington-Finkelstein and Painleve-Gullstrand coordinates 

For the present article we will use geometrical units such that the speed of light c and 
Newton’s gravitational constant G are set to unity. On top of that, it is convenient 
for the following analysis to measure all radial distances in units of the Schwarzschild 
radius r s = 2M. With this convention the horizon of the Schwarzschild black hole is 
located at r = 1 and the light-sphere, where circular lightlike orbits exist, is found at 
r = 3/2. The Schwarzschild metric, written in the usual Schwarzschild coordinates, 
becomes then 

gu V dx^dx v = — — — —dt 2 + - —-— -dr 2 + r 2 (dd 2 + sin 2 ddp 2 ) . (1) 

r (r — 1) v ' 

As an alternative to the usual Schwarzschild coordinates (t , r, d, ip) we can introduce 
other coordinate systems to describe the same spacetime geometry. Two examples 
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that we employ in this work are the Eddington-Finkelstein (EF) and the Painleve- 
Gullstrand (PG) coordinates which are both regular at the horizon. For details on 
regular coordinate systems for the Schwarzschild spacetime we refer the reader to 
Ref. |3]and references therein. The ingoing version of the EF coordinate^ is obtained 
using 

vef = t + r + log(r - 1), (2a) 

u = r , (2b) 


and the Eddington-Finkelstein coordinate v is constant along all radially ingoing 
lightlikc geodesics. The Painleve-Gullstrand coordinates El are in our units given by 


vpg = t + 2y/r + log 


\fr - 1 


y/r + 1 


u = r. 


(3a) 

(3b) 


and are associated with infalling observers that follow radial timelike geodesics 
starting from spatial infinity P Both coordinate systems are obtained by a coordinate 
transformation starting with usual Schwarzschild coordinates t and r, while the 
angular coordinates are not changed to preserve the symmetry of the situation. The 
Schwarzschild metric is not diagonal anymore in these coordinate systems but the 
apparent singularity at the horizon can be characterized as being only a coordinate 
singularity since the metric is invertible at r = 1 in EF and PG coordinates. In the 
following, we will omit the subscripts “ EF ” and “PG” at the new time coordinate 
v since the two distinct cases are treated in separate sections. 


2.2. The confluent Heun equation 


The singly confluent Heun equation (CHE), which belongs to Heun’s class of differ¬ 
ential equations,®® can be given in the canonical form 


_d^_ 

dz 2 


6+1 c + 1 


z — 1/ dz 


- 7 -+ - 


z-1 


y(z) = o, 


( 4 ) 


where the two parameters d and e are defined in terms of all the other ones by 


/i = — (a — b — c — 2e + a6 — 6c) , 
v = - (a + b + c + 2d + 2e + ac + be) . 


(5a) 

(5b) 


a In 1924 Eddingtoi® was the first to discover this coordinate system which demontrates that the 
singularity at the horizon is a coordinate artifact and not a physical singularity, even though he 
never stated this explicitly. The same coodinate system was rediscovered by D. Finkelstein in 
1958® 

b P. Painleve (1921) and A. Gullstrand (1922) independently discovered a spherically symmetric 
metric without realizing that it was just the Schwarzschild metric in new coordinates. 
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The form is also called the Maple-form of the CHE since it is implemented 
in the modern computer algebra system MapleP The CHE ® possesses regular 
singularities at z = 0,1 and a single irregular singular point at z = ooP® 

At the regular singularities local solutions can be obtained in terms of Frobenius 
series 10 or be given using the standard confluent Heun function HeunC that is im¬ 
plemented in Maple as well. For details about the definition of this confluent Heun 
function the reader is referred to, e.g., Refs.[l]j9] At the irregular singularity asymp¬ 
totic series, also called Thome-type solutions, can be constructed! 8 The behavior of 
the local solutions depends on the characteristic exponents of the considered sin¬ 
gularity. For details on the construction of these local solutions, that involves the 
so-called generalized Riemann scheme (GRS) as a very useful tool, the reader is 
referred to Ref. Q] and the excellent book on singular differential equations [8] 


3. Massive scalar field in EF coordinates 


In this section we complement the recent work in Ref. (T] and give the solutions of 
the perturbation equation for a massive scalar field on Schwarzschild background 
spacetime in EF coordinates. The propagation of such a massive scalar field is 
described by the Klein-Gordon equation with a mass term 


(□ — TO 2 ) $ = — to 2 ) $ 


dx» 






where to is the mass of the scalar boson. The scalar Klein-Gordon field $ is a 
function of all the coordinates t, r, i?, ip and can be expanded into scalar spherical 
harmonics Yi m ('d ) ip) as shown in many standard textbooks, see e.g. Ref. [Til 


oo L -p /. \ 

$ = $,¥>)= ^ ^2 ———Y lm {'d,p) , 


( 7 ) 


Z=0 m=—l 


where we introduce the tortoise coordinate 


r* = r + log(r — 1). 


( 8 ) 


With the ansatz we obtain the wave equation 


drl 


dt 2 


Vi(r) 


^i(t,r *) = 0, 


Vi(r) 


(r- 1) 


( 1(1 + 1 ) 

V r 2 



(9) 


This equation is sometimes called the time-dependent Regge-Wheeler equation for 
a massive spin s = 0 field. In their original work Regge and Wheeler^ derived the 
analogous equation for the massless s = 2 case. The lowest mode (l = s = 0) and 
the first excited mode Z = s + l = lof the potential Vj(r) are shown in Fig. 0 for 
different values of the mass to. The maxima are located near the light-sphere radius 
r = 3/2. In the limit r — > oo the value of the potential is given by the squared mass. 
This property distinguishes the present situation from the massless case which is 
obtained in the limit to —> 0. A detailed numerical analysis of (quasi-) resonant 
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Fig. 1. The potential as a function of the tortoise coordinate r**: the lowest mode l = s = 0 
(left) and the first excited mode l = 1 (right) for m = 0.25 (dotted), m = 0.3 (dashed), m = 0.35 
(dash-dotted) and m = 0.4 (solid). 


states (i.e., the gray region shown in Fig. [T]) in this potential can be found in Ref. 
EJ where Schwarzschild coordinates were used in the discussion but EF coordinates 
were considered for the numerical evolution as well. 

Using the results of Ref. CD we can transform the wave equation ([9]) into a wave 
equation in Eddington-Finkelstein coordinates that allows us to obtain analytical 
and regular solutions around the black hole horizon at r = 1. There is no known 
exact analytic solution for the wave equation in the form ([9]) due to the involved 
potential term. The corresponding wave equation in EF coordinates iJ2 


i — 1 


dr 2 


(r - 1) d 2 (r — 1) d 


drdv 


dr 


~ Vi(r) 


*i{v,r) = 0. (10) 


This equation can be further transformed into a radial equation using the ansatz 

*(«,r) =e- iuv R ul (r) (11) 


that separates the EF time coordinate v. The radial equation that the function 
R u i(r) has to fulfill is then given by 


[ d2 +1 

(1 — 2 icor 2 

d i 

[dr 2 

V r(r - 1) ) 

dr V 


+ 


1 


+ 


m 2 r 


( Kl + 1) 

\r(r — 1) ' r 2 (r — 1) ' (r - l) J_ 


Rci(r)= 0. (12) 


We recognize regular singularities at the origin r = 0 and at the horizon r = 
1. Moreover, there is an irregular singularity at z = oo which can be seen by 
substituting ( = 1/r and investigating the behavior at ( = 0. With these properties 
the radial equation belongs to the class of singly confluent Heun equations. The 
transformation 


R ujl ( r ) = re( iu + Vm ^>H ull (r) 


( 13 ) 
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yields a CHE in the form of dU) with the parameters 

a = 2\J m 2 — w 2 , 6 = 0, c=— 2iui, fj, = iui + 1(1 + 1) + \/ to 2 — co 2 , 

v = —2 i\Jm 2 — uj 2 uj — l(l + 1) — m 2 + \/m 2 — oj 2 + 2 w 2 — ico . (14) 

Hence, the GRS for the radial equation becomes, using the general result of Ref. _lj 


(l 1 2 

0 1 oo 

0 0 (1-2 iuj + 


2 w 2 — m 2 


0 2 ico 


1 - 


2y/m 2 — w 2 
2w 2 — m 2 


V 


(15) 


/ 


2 \fm 2 — oj 2 , 

0 

(—2 \/m 2 — w 2 ) 

We can give two independent solutions i?^) /7 (r; 1) at the regular singularity at the 
horizon (r = 1) in terms of the standard confluent Heun function 

Rh{ r -l) =r 

x HeunC(—2 \Jm 2 — w 2 , 0 ,to 2 — 2 + 1) + 2w 2 — m 2 ,1 — r), 

(16a) 

R^(r; 1) = r e( iaj+Vm "-^) r (r - l) 2itJ 


x HeunC(—2 \Jm 2 — w 2 ,2iw, 0, m 2 — 2w 2 , — Z(Z + 1) + 2w 2 — m 2 ,1 — r). 

(16b) 


From this we recover indeed the previous solutions of the massless casd^ in the limit 
to —> 0. While the first solution 1) is regular at the black hole horizon and 

corresponds to purely ingoing radiation at that radius, the second solution l?^((r; 1) 
is not regular there since the phase performs infinitely many turns on the unit circle 
in the complex plane due to the term (r — 1) 2 *“ that has no well defined limit for 
r —> 1. When we want to imply causal boundary conditions, i.e., no radiation should 
escape the black hole horizon, we only need to consider the first, regular, solution. 
Figure ([2]) shows the l = 0 to l = 5 modes of this regular solution 1). 

The local solutions (1161) are, unfortunately, only useful in the rather small domain 
of convergence |r — 1| < 1 but can be analytically continued to the entire domain 
r £ ]0, oo[ by the procedure that was demonstrated in Ref.Q]for the massless case and 
a general bosonic spin s in EF coordinates. We will sketch how such a continuation 
procedure works for the case of massless fields in PG coordinates in the next section. 


4. Massless fields in Painleve-Gullstrand coordinates 

We start now with a brief reconsideration of the results derived in Ref. [1] that will 
be needed for this section. In that article we have shown that the wave equation on 





















May 28, 2015 0:31 WSPC/INSTRUCTION FILE PhilippPerlickS'ijmpd 


Schwarzschild radial perturbations in EF and PG coordinates 7 






r 




r 


Fig. 2. The regular solution 1) for a massive scalar field in EF coordinates for Z = 0,1 

(top), l = 2,3 (middle), / = 4, 5 (bottom) for the choice to = 4, m = 0.4. The modulus (solid) and 
the real part (dashed) are shown. The horizon of the Schwarzschild black hole is located at r = 1. 


Schwarzschild background that describes the propagation of massless scalar (s = 
0), electromagnetic (s = 1) or gravitational (s = 2) perturbations written in PG 
coordinates is given by the partial differential equation 



2 


dr 2 


(r — 1) d 2 9 (r — 1) d 2 
r dv 2 r 3 / 2 drdv 

(r — 1)9 (r — 1) d 
2r 5 / 2 r 3 


V.i(r) 


^(v, r) = 0 , 


(17) 


where the spin-dependent potential is given by 


V s i(r) 


(r ~ 1 ) 

r 3 


1(1 + 1 ) 


1 - s z 


r 


(18) 
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We can separate now the PG time coordinate v using the separation ansatz 

V(v,r)=e- i “ v R usl (r). (19) 

Thereby we introduce the frequency w of the perturbation that has the dimension 
of an inverse length in our units; so the combination loM is a dimensionless quan¬ 
tity. The harmonic time dependence in the ansatz © leads to a purely radial 
perturbation equation, that is 




+ I r 2 uj 2 + 


luir 


1/2 


- 1 ( 1 + 1 )- 


1 — s* 


R U si(r) = 0. ( 20 ) 


This radial equation is a second order, ordinary, linear and homogeneous differential 
equation. We further recognize that equation (EU1) possesses regular singularities 
at r = 0,1. Using the transformation £ = 1/r we are able to find an irregular 
singularity at £ = 0, which correspond to r = oo. With these singularities the radial 
equation could belong to the class of singly confluent Heun equationsP® Indeed, a 
transformation that maps (l20l) into the Maple-form of the confluent Heun equation 
is given by 


R us i(r) = r s+1 (Vt-ir(r 3 ' 2 -Vt + r-l) *“ e “< 2 ^ +r >ff us! (r). (21) 

The function H^si (?') has to fulfill the differential equation (U) with the parameters 

a = 2 ioj , b = 2s , c = —2 ico , /r = 4 iujs + 2 iu + 1(1 + 1) — s(s + 1), 
v — — 2/iujs + 4 uj 2 + s(s + 1) — 1(1 + 1) ■ (^^) 

Hence, the GRS for this equation can be given by using the general result,^ which 
yields 

/ 1 1 2 \ 

0 1 oo ; r 

0 0 s + 1 — 2 iu 

—2s 2 iu s + 1 

0 

y —2 ico 


(23) 
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Fig. 3. The two local solutions around the horizon: the regular solution j(r; 1) (left) and the 
second solution R J J s i (r; 1) (right) for uj = 6, l = s = 0. The modulus (solid) and real part (dashed) 
are shown and the horizon is located at r = 1. 


4.1. Local solutions at the horizon 

We can construct two independent solutions at the regular singularity at the horizon 
(r = 1) using the confluent Heun function.^® These two solutions are 


( , \ —ito . 

r 3/2 — \Jr + r — l) e M 2 v^+r) 

x HeunC(— 2iu, — 2iw, 2s, —2 to 2 , —1(1 + 1) + s 2 + 2 lo 2 , 1 — r), (24a) 

R T J sl (r; 1) = r s+1 (v^ - l) iw (r 3/2 - VF + r - l) e ” (2 ^ +r) 

x (r — l) 2l “HeunC(—2zw, 2*w, 2s, —2w 2 , —Z(Z + 1) + s 2 + 2w 2 ,1 — r). 

(24b) 


The first solution R^ s i(r', 1) describes a wave that propagates into the black hole 
horizon and is totally regular when approaching this radius. The second solution 
R/Jsi (r; 1) on the other hand describes outgoing radiation at r = 1 and is not well 
defined when approaching the horizon since the phase of this solution performs 
infinitely many turns on the unit circle in the complex plane in the limit r —> 1, as 
we have seen in section [3] before. Figure <[3|) shows the real part and the modulus of 
both solutions and clearly demonstrates the corresponding behavior at the horizon. 


4.2. Asymptotic series at spatial infinity 

Using the entries of the third column of the GRS (l23l) we can construct the Thome- 
type solutions (asymptotic series) at the irregular singularity at spatial infinity. 














May 28, 2015 0:31 WSPC/INSTRUCTION FILE PhilippPerlickS'ijmpd 


10 Dennis Philipp and Volker Perlick 
These solutions ar(P® 

RLl (r; oo) = r 5+1 (VF - 1)“ (r 3 / 2 -^r + r- 

OO 

^ g iuj(2y/r-\-r ) E p k r~ {k+s+1 ~ 2iul) , 

k =0 

Rliii^o o) = r s+1 (v/f- 1)“ (r 3/2 - Vr + r- 

OO 

x e i^(2^+r-) e -2io;r ^ 0 - fcr --( fc + s + 1 ) , 
fc=0 

and the series expansion coefficients p*, and tife can be calculated in terms of recur¬ 
rence relations that are obtained by inserting the solutions above into the differential 
equation. The first solution, together with the separated time dependence ~ e ~ tulv 
corresponds to radiation that is propagating into the singularity at r = oo (outgo¬ 
ing), while the second solution describes waves that emerge from this singularity 
and propagate to smaller radii (ingoing). 

4.3. Analytic continuation of Frobenius solutions 

Both local Frobenius solutions at the horizon R : pJ (r; 1) are convergent within the 
unit circle centered at r = 1 in the complex plane as explained in Ref. [T] However, 
it is possible to perform an analytic continuation that combines the Frobenius so¬ 
lutions (12411 with the asymptotic series (l25l) and in this way to obtain solutions on 
the entire open interval r €]0,oo[ . This procedure can be partly traced back to 
work of Jaffe 13 and is demonstrated in detail for the radial perturbation equation 
for massless fields in EF coordinates in Ref. H] using the general formalism intro¬ 
duced by Slavyanov and Lay 8 in terms of the central two point connection problem 
(CTCP). We refer the reader to these works and restrict ourselves here to a brief 
sketch of how to apply this procedure to the present case of PG coordinates. The 
five-step procedure that we apply is the following: 


(25a) 


(25b) 


1. Shift of singularities and reordering the differential equation 


We perform the simple transformation z = r — 1 to shift the involved regular 
singularity from r = 1 to z = 0 and reorder the resulting differential equation in 
the form 


’il 

dz 2 


■ 2* 




C 2 

z + z* 



yu,3i(z) = 0 


(26) 


with the parameter 


z* = 1, Go = 2iu , A\ = 1 — 2 ioj , Ai = 1 + 2s, 

C\ = — 2iajs -(- 4w 2 -}- s(s + 1) — l{l -\~ 1), C 2 — diojs 2iuj T l{l ffi 1) — s(s -(- 1), 

(27) 
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to match the general formulas given in Ref. [8] Note that the irregular singularity 
at r = oo remains at z = oo. 


2. s-homotopic transformation of the dependent variable 

The next step is to perform a special transformation of the dependent variable 
H us i(z). The transformation we need is a so-called s-homotopic transformation and 
given by 


UqjsI^Z) I y ZZcjsliz) 

H usi {z) = e vz z^ (z + 1 r E U 3 i(z) . (28) 

The three new exponents iz, p i,P 2 are connected to characteristic exponents at the 
two involved singularities. Similar to the treatment in Ref. [1] we obtain the possible 
values 


- pi is the indicial exponent of the chosen Frobenius solution at z = 0. 
=» pi = 0 or pi = 2 iuj according to the GRS K51) . 

- v is the characteristic Thome exponent of order one and depends on the 
chosen Thome solution at z = 00 . The possibilities are v = 0 or v = —2 iui. 

- For large 2 we get {z+ 1) M2 ~ Thus, pi+p -2 is the characteristic 

exponent of order zero of the chosen Thome solution. => p 2 = (s + 1), 
P 2 = (s + 1) — 2 ioo or p 2 = (s + 1) — 4 iw. 

As a summary of this step we sketch in Figure d4j) the appropriate choice of the 
three exponents according to the behavior at the singularities that is determined 
by the Frobenius solutions and the asymptotic series. 


Frobenius solutions 


Ks,(r; 1) Pl = 0r~ 


/12 = (s + 1) - 2 iu 


Thome solutions 
— v = ~ 2iu] R Li(r;oo) 


X.V. 


X 


y 

a:/ 

❖V 


Rl Jsi(r; 1 ) 


Pi = (a + 1) - 2 iuj 


r = 0 Rfiii r;oo) 


Fig. 4. A scheme showing the proper choice of the exponents is, and 112 related to local solutions 
for the connection problem. 
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3. Mobius transformation of the independent variable 

We apply now a transformation of the independent variable. The necessary trans¬ 
formation in this step is given by the Mobius transformation 

241= ——— (29) 

2 + 1 

that maps the irregular singularity from z = oo to x = 1 at the boundary of the unit 
circle. The other singularity involved is now located at x = 0 and we have mapped 
the infinite interval 2 € [ 0 , oo] to the finite interval of unit length x € [ 0 , 1 ]. 


4- Series expansion ansatz 

Now, we solve the resulting differential equation for {x) with the series expan¬ 
sion ansatz 

OO 

^ S [(i) = ^(t(«,s,l) i 1 . (30) 

fc=o 

This yields a three term recurrence relation for the series expansion coefficients 
£fc(w, s, /) in the form 

Co = arbitrary 

«oCi + A)£o = 0 

QJfcCfe+l + Pkfk + 7fc£fc-l = 0 , (31) 

where the a*,, /3k and 7 k depend on the parameters (1271) and, therefore, also on w, s 
and l. The actual form of these quantities can be derived using the results given in 
Ref. |T] The series representation of the function H ws ; (x) is absolutely convergent on 
the interval |x| < 1.EE3 


5. Full solution of the radial equation 


To construct the total solution to the radial equation in PG coordinates (l20l) we 
have to recover the radius variable r that we started with. This is done by using 
the inverse Mobius transformation 


x = —— = — . (32) 

2 + 1 r 

Now, we construct the total solution using again the transformation Hence, 
we obtain 


( , \ —ioj , _ 

r 3 / 2 — \fr + r — lj e M^+r) 

OO / -j \ fc 

x e^-^r - lfr^ f fc (w, a, , (33) 

fc =0 

where the exponents pi t 2 and v must be chosen according to the physical bound¬ 
ary conditions with the help of Figure (j4| and the follow from the mentioned 
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three-term recurrence relation. The solution (1331) is now absolutely convergent in 
the interval r s]0.5, ooP and we have, thus, solutions on the the entire interval 
r e]0,oo[ since for the region r e]0,2[ the Frobenius solutions (l24l) in terms of 
the standard confluent Heun function can be used and the analytically continued 
solutions coincide with these local Frobenius solutions in the overlapping region. 


5. Conclusion 

We have shown that the propagation equations for bosonic spin s fields on 
Schwarzschild background spacetime can be reduced to purely radial differential 
equations for Eddington-Finkelstein and Painleve-Gullstrand coordinates. In these 
coordinate systems the radial equations admit (local) solutions in terms of conver¬ 
gent series. We complemented our recent workf^ where massless fields were inves¬ 
tigated in EF coordinates, by additionally analyzing a massive scalar field in EF 
coordinates and extending the treatment of the massless bosonic fields to the case 
of PG coordinates. In each of these cases it is possible to derive two local solutions 
around the regular singularity at the black hole horizon. One of these solutions is 
regular at r = 2 M and connected to causal boundary conditions of purely ingo¬ 
ing radiation, and the other solution describes radiation emerging from the black 
hole. While the regular solution can be used to investigate scattering of radiation or 
to calculate quasi-normal modes, the second solution describes Hawking radiation 
and is, therefore, connected to black hole evaporation as shown in Ref. 21 We have 
further shown that using an analytic continuation procedure the bounded domain 
of convergence of the local solution (i.e., the unit circle around the horizon in the 
complex plane) can be expanded to the entire open interval between the origin and 
spatial infinity. It is a remarkable feature that the structure of the radial equations 
in EF and PG coordinates is the same. Both cases belong to the class of confluent 
Heun equations and can be solved using Frobenius series at the regular singularity 
at r = 2 M that can analytically be continued up to, but not in general including, 
spatial infinity. The equations in PG coordinates are somewhat more complicated 
and we recommend, thus, to treat all these perturbations in EF coordinates. If, 
however, the situation should be linked to the measurement of a radially infalling 
observer on a timelike geodesic, one should use the solution we provided in Eqn. 
m■ This solution can be used to investigate various phenomena concerning mass¬ 
less bosonic waves on Schwarzschild background and can be adopted to a broad 
range of initial or boundary value conditions. The most interesting feature of the 
regular coordinate systems is the possibility to penetrate the horizon and thereby 
to follow the radiation that falls into the black hole. 
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